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ABSTRACT 

The rare phenomenon of ribbon-like hard X-ray (HXR) sources up to 100 keV found in the 2005 May 
13 M8.0 flare observed with the Reuven Ramaty High Energy Solar Spectroscopic Imager provides detailed 
information on the spatial distribution of flare HXR emission. In this Letter, we further investigate the charac- 
teristics of HXR emission in this event using imaging spectroscopy, from which we obtain spatially resolved 
HXR spectral maps during the flare impulsive phase. As a result we found, along a flare ribbon, an anticorre- 
lation relationship between the local HXR flux and the local HXR spectral index. We suggest that this can be 
regarded as a spatial analog of the well-known temporal soft-hard- soft spectral evolution pattern of the inte- 
grated HXR flux. We also found an anticorrelation between HXR spectral index and local electric field along 
the ribbon, which suggests the electron acceleration by the electric field during flares. 
Subject headings: Sun: flares — Sun: X-rays, gamma rays 



Oh; 

o : 
a: 



> 

o 

m 



o 



1. INTRODUCTION 

It has long been recognized that hard X-ray (HXR) emis- 
sion is a powerful diagnostic of accelerated energetic elec- 
trons produced by flares. The HXR spectrum emitted by non- 
thermal electrons often appears as a power-law distribution in 
photon energy, which implies a characteristic in the energy 
distribution of the electron flux bombarding the target under 
the bremsstrahlung emission mechanism. Temporal evolu- 
tion of the electron energy distribution is highly important for 
identifying the dominant acceleration process. 

As known from the early results of the traditional 
scintillation-counter spectrometers, the HXR spectral index 
generally follows a soft-hard- soft (SHS) spectral pattern in 
the rise-maximum -decay phase of flares ([Parks & Wi nckled 
[T^:lBen3fT977l: iBrown & Loranl fl985l: lDennislll985lT lt 
was further corroborated by the results of Reuve n Ramaty 
High Energy Solar Spectroscopic Imager (RHESS 7: lLin et all 
2002) tha t can measure HXR spe ctra with much higher res- 
olution. iHudson & Far nik ( 2002) illustrated the consistency 
of the SHS pattern derived using higher energy resolution 
RHESSI data (of order 1 keV) with that obtained with the 
HXRS scintillation-co unter spectra. A systematic study by 
iGrigis & Benzl (120041) using RHESSI revealed that the SHS 
pattern appears even in sub-peaks of HXR with durations of 
one minute to shorter than 8 s. With RHESSFs capability of 
imaging spectroscopy, the study of SHS has been exte nded 
to spatially resolved HX R sources. Ems lie etail (l2Q03h and 
Battaglia & Benz ( 2006) examined spectral evolution of both 
coronal and footpoint HXR sources in the course of flares and 
find that they commonly show the SHS behavior. Based on 
the results they suggested that SHS can be an intrinsic feature 
of the electron acceleration in the flare impulsive phase. 

The SHS can also be regarded as a relation between HXR 
spectral index and flux because the ha rdest spectrum usuall y 
appears in the period of maximum flux dGrigis & Benz 2006). 
As a possible extension of this relation to spatial characteris- 
tics, we may consider a spatial analog of the SHS pattern in 
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which a stronger HXR region has a harder spectrum. Such 
idea has been propo sed in a couple of studies. For instance, 
Mas uda et alJ (l200l|) determined spectral index maps of an 
extended HXR event using Yohkoh data to find spectral hard- 
ening toward the rib bon edge where the m ost intense energy 
release is expected. IHudson et al.l (|2004) proposed that the 
regions of weaker HXR emission would correspond to softer 
HXR spectra, in an attempt to explain the footpoint-like HXR 
morphology in contrast with extensive ribbon structure seen at 
Ha and UV wavelengths. However, it will need a more thor- 
ough examination to see whether the relation between HXR 
flux and spectral index in time can simply be transformed to 
such a spatial analog. This is one of the goals of this paper that 
we intend to achieve using RHESSI imaging spectroscopy. 

In addition to spectral evolution, the HXR flare ribbon mo- 
tion can also manifest the progression of flare energy re- 
lease due to magnetic reconnection. This was enabled after 
iForbes & Priestl (1 19841) derived a relationship between the rib- 
bon expansion velocity u and the electric field E in the recon- 
necting current sheet (RCS) in the form of E = uB, where B 
is the local vertical magnetic field strength in the footpoint. 
It is important to have an indirect measurement of the elec- 
tric field, because direct acceleration of electrons by the elec- 
tric field is a candi date mechanism f or creating high energy 
electrons in flares (lLitvinenkol fT996). Some studies showed 
that there is a temporal correlation betw een HXR flux and 
the electric field derived this way (e.g.. lOiu et al.ll20()2|) . It 
is then of next interest whether there is not only a temporal 
but also a spatial correlation betwe en HXR spectral index and 
electric field in the RCS. Recently. Wood & Neukirchl (120051) 
presented a particle simulation with prescribed electric and 
magnetic fields, which shows a hardening of the electron en- 
ergy spectrum with increasing electric field. Such a model 
prediction can be compared with the above mentioned mea- 
surements of the HXR spectrum and electric field. 

In this Letter, we investigate spatial distribution of HXR 
spectrum and electric field duri ng the 2005 Ma y 13 M8.0 
flare. This is a particular event dLiu et al.ll2007b|) . where the 
HXR sources are so extended as to be suitable for the imag- 
ing spectroscopy, and is thus adequate for studying the possi- 
ble spatial relationship between the spectral index and electric 
field. 
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Fig. 1. — RHESSI CLEAN images for five one-minute time intervals in 
four energy bands from 9-100 keV. The images are centered on the point 
(-168", 230") with solar west to the right and north up. The field of view 
is 128" x 128". The white dashed line superposed on the lower left panel 
denotes the magnetic polarity inversion line. 



2. IMAGING SPECTROSCOPY 
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RHESSI imaging spectroscopy has thus 
in such a way to determine the flux 
over spatially isolated sources (e.g , 
Battaglia & Benz 2006). This technique gives a viable result 
when the source is so compact as to be barely resolved. In 
our case, we however study extended HXR sources and our 
goal is to explore the spatial variation of spectral index inside 
the sources. We thus choose to measure local HXR spectral 
index based on spectrum in each pixel. Previously, such at- 
tempt was also made for the only o t her r eported ribbon-like 
HXR event, in which iMasuda et al.l (120011) obtained spectral 
maps providing the spectral index information in each pixel, 
from the count ratio between Yohkoh M2- and H-band. With 
RHESSFs higher temporal and energy resolution as well as 
improved dynamic range, we expect a better imaging spec- 
troscopic result. For this moderate event, however, we still 
suffer from the trade off between choosing finer energy/time 
bins and making better images by accumulating for wider en- 
ergy/time bins. 

Most often, RHESSI imaging is mad e using the CLEAN 
algorithm with grids 3-9. Recently, iDennis etaD |2007) 
showed that adding grids 1 and 2 (the finest grids) can enhance 
the image quality depending on the source structure, and the 
ribbon-like HXR sources in this event is such a suitable case. 
We found that use of all grids not only improves the imaging 
quality compared with using grids 3-9, but also gives bet- 
ter count statistics that is essential for imaging spectroscopy. 
We also found that natural weighting produces a m ore phys- 
ically r easonable source morphology for this event dLiu et al.l 
I2007ah . While PIXON or NJIT-MEM may give more details 
in the source morphology, CLEAN gives a reason able result 
for imaging spectroscopy dBattaglia & Benzl [20061) . We thus 
use the CLEAN algorithm along with all grids and natural 
weighting, which gives a FWHM resolution of ~5.9". 

We first make 128 x 128 pixels images with a pixel size of 
\" for the five one-minute interva ls that cover the rise and 
decay phases of HXR (a-e as in Liu etal.ll2007ai Fig. 1). 
Figure [T] shows the images in four energy bands from 9- 
100 keV. The emission at high energies (>25 keV) are mainly 
sources located in each side of the magnetic polarity inver- 




FiG. 2. — HXR spectral maps in the five one-minute time intervals in Fig.[T] 
covering the HXR rise to decay phases. These maps are generated based 
on HXR images in Fig. [T] assuming a single power-law spectrum at each 
pixel. A lower power-law index corresponds to a flatter spectrum. The white 
background indicates the position where no good power-law fitting can be 
obtained. The black curves are our trace-outs of multiple locations along Ha 
ribbons determined by finding out local maximum Ha intensity in each time 
interval. The curves are composed of 100 points evenly distributed along the 
ribbon. 



sion line (PIL) presumably from the footpoints of the flaring 
loops. A footpoint-to-ribbon transformation of the HXR mor- 
phology can be clearly seen across the peak of the HXRs at 
16:42:04 UT. At lower energies, sources above the PIL be- 
come evident and they could come from the tops of the loops 
joining the HXR footpoints and ribbons, while the footpoint 
emissions are still visible, as low as in 9-12 keV band. We 
emphasize that the ribbon-like HXR sources are an intrinsic 
feature in this ev ent, as evidenced by its fidelity to t he UV 
ribbon emissions dLiu et al . 112007 al:lDennis et al.ll2007l) . 

We then derive spectral index maps from the four energy 
bands shown in Figure [TJ assuming that the local HXR spec- 
trum in each pixel follows a power-law distribution in pho- 
ton energy with spectral index defined by 7 = -dlogl/dloge. 
While it is more desirable to have narrower energy bands for 
a better spectroscopy, the choice of only four energy band is 
inevitable to ensure enough photon statistics in each time in- 
terval for this event. Another issue in HXR imaging spec- 
troscopy is that the coronal and footpoint sources show dif- 
ferent spectral characteristics. Coronal HXR sources usually 
show ther mal-like spectrum and f ootpoint sources, power-law 
spectrum (Battagli a & Benzl 120061) . and we had better distin- 
guish them from each other. In the present event, the HXR 
emission is dominated by the footpoint sources in a wide en- 
ergy range ~10-100 keV and we fit all sources in the energy 
range to power-law spectra. Nevertheless, our approxima- 
tion that each photon spectrum measured in four energy bands 
should follow a single power law may not be good enough in 
all pixels. We thus set a threshold, based on the chi- square 
of a value of 0.25 returned by the LINFIT procedure of IDL, 
for the goodness of the fitting, below which we do not cal- 
culate the spectral index. The results for each time interval 
are shown in Figure [3 where the color coded spectral index 
is displayed only for those fittable pixels. To ascertain that 
the flux in each pixel is enough for generating a reasonable 
spectrum, we further reduce the images to 32 x 32 pixels by 
averaging a set of 4 x 4 pixels into one. We repeat the fitting 
procedure and obtain nearly identical spectral features. 

3. RESULTS 
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Fig. 3. — Spatial distribution of HXR spectral index (black) in comparison with those of HXR flux in the 25-100 keV range (grey, upper panels) and electric 
field (grey, lower panels) in the five one-minute time intervals, shown as functions of the ribbon distance. Depicted as the black line in Fig. [2] the index of the 
ribbon distance runs from in the northern end to 99 in the southern end. The existence of break points in the profiles of spectral index indicates bad fittings at 
those locations. 

Comparison of Figure [T] and Figure [2 reveals several char- 
acteristics. First, there are three major footpoint-like HXR 
emitting sources with hard spectra near the HXR peak at 
^16:42:04 UT, when the spectral index reaches a mini- 
mum value of ~2.2. The averaged value of spectral index 
over the fittable pixels in the field of view is ~3.3 at this 
time, which agrees with that derived using the OSPEX pack- 
age for integrated X-ray emission. Earlier and later in the 
event, the overall spectrum of the flare region as well as 
those of the main HXR sources are seen to be much softer, 
with higher index values. This temporal evolution of the 
SHS pattern is thus wh at typically has been observed before 
(Battaglia & Benz 2006). Second, the main HXR sources 
show a spatial distribution of spectral index from the cen- 
ter with smaller value (harder spectrum) to the outer regions 
with larger value (steeper spectrum), which is most promi- 
nent near the HXR peak (see panel 16:41:34-16:42:34 UT). 
When the HXR evolve to a ribbon morphology later in the 
event (e.g., panel 16:43:34-16:44:34 UT), there are still ker- 
nels with harder spectra discernible, although the whole sys- 
tem has a much steeper spectrum compared with the flare 
peak. 

As our major interest lies in how the physical quantities 
vary spatially, we make the comparison between the HXR 
spectral index and flux along the ribbon axis. Specifically, we 
measure the variation of spectral index along the eastern flare 
ribbon, where there is a clear footpoint- to-ribbon evoluti on of 
HXR morphology (c.f. Fig.[TJ also see lLiu et aLll20 07a). We 
trace out the spectral index and the flux in the 25-100 keV 
ra nge using the sam e indexes of ribbon distance as defined 
by I Jing et all (I2007L specified in Fig. 2 as black line at each 
time interval), and present the results in the upper panels of 
Figure [3j It is obvious in each time interval that the HXR 
spectral index exhibits a strong spatial anticorrelation with the 
HXR flux (log/ 7 ), with absolute values of the correlation co- 
efficient >0.8. In the lower panels of Figure[3l we compare 
the spatial evolution of the HXR spectral index with that of 
the electric field (logE) in the RCS, which was previously de- 
rived by tracing the Ha ribbon motion an d incorpor ating the 
longitudinal magnetic field measurement ( Jing et al. 20Q7|). It 
can be seen that there also exists a prominent anticorrelation 
relationship between these two quantities, with absolute val- 
ues of the correlation coefficient >0.65. 

Figure [4] shows the results in Figure [3] as scatter plots. The 
relationship between the spectral index and flux and the spec- 
tral index and electric field are shown in the upper and lower 
panel, respectively. First, we can see without ambiguity that at 
a specific time, the source position with a weaker HXR emis- 
sion (lower HXR flux) correspond to steeper X-ray spectra 
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Fig. 4. — Scatter plot of spectral index (7) and HXR flux (log/ 7 , upper 
panel), and spectral index and electric field (log£, lower panel), measured at 
each indexed location along the ribbon for all the time intervals. Data points 
and the best fits in linear-log space (solid lines) by minimizing the chi-square 
error statistic for different time intervals are represented with different colors. 



and presumably to softer electron precipitation spectra. We 
therefore suggest that this anticorrelation between HXR spec- 
tral index and flux is a spatial analog of the well-known tem- 
poral SHS pattern of HXR emission. Second, the anticorrela- 
tion relationship between the spectral index and electric field 
clearly points to a softer HXR spectrum in the case of weaker 
electric fiel d, consistent with the tren d predicted by numerical 
simulation (I Wood & Neukirchl2005|) . This result thus appears 
to support the hypothesis that direct acceleration by the elec- 
tric field in the RCS may play an important role in producing 
energetic electrons in flares. 

4. SUMMARY AND DISCUSSION 

Exploiting RHESSFs capability of imaging spectroscopy, 
we have determined spectral index of the local photon spec- 
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trum in the unusually extended HXR source observed in the 
2005 May 13 M8.0 flare. We then find a spatial anticorrela- 
tion relationship between both the local HXR flux and electric 
field corresponding to the local HXR spectral index at several 
time intervals. We discuss the present results in comparison 
with similar works on solar HXR imaging spec troscopy. 

The present approach is closest to that of iMasuda et al.l 
( 2001) who used Yohkoh images at two energy bands to find a 
spectral index change across flare ribbons. In their result, the 
HXR spectrum at the outer edge of the ribbon is found to be 
harder than in other regions of the ribbon and the hardness of 
the spectrum gradually changes across the ribbon width. This 
is consistent with the physical picture that the ribbon edge 
is connected to the most recently reconnected field lines and 
thus shows the most energetic electrons. The present study 
shows a similar result with more spatial details. It clearly 
shows the hardest HXR sources lying a long th e edge of flare 
ribbon in UV (c.f. Fig l2l and IQu et al.ll2007aL Fig. 2) and a 
smooth transition of spectral index across the HXR sources. 
Particularly, we find that the HXR spectral index exhibits a 
strong spatial anticorrelation with the HXR flux in all time in- 
tervals during the flare impulsive phase. We call this spatial 
SHS behavior in analogy with the well-known temporal SHS 
pattern of integrated HXR flux. The spatial SHS may also 
be an essential feature of solar flare electron acceleration, and 
may help explaining the confined nature of HX R sources com- 
pared with extended Ha and UV ribbons (c.f. Hudso net al.l 
12001 . 

We must note that the spatial SHS implies a more strict rela- 
tionship between HXR flux and spectral index than that found 
in the previous RHESSI stud ies on the temporal SHS behavior 
in ind ividual sources (lEmslie et al.l 120031: iBattaglia & Bend 
2006). In the latter results, each isolated source exhibits the 
temporal SHS pattern independent of each other and thus the 
normalization of the flux-spectral index relation (F-S rela- 
tion) may differ from one source to another. On the other 
hand, our spatial SHS implies that the same normalization 
of the F-S relation should hold in all regions. A major dif- 
ference between those works and ours lies in that our spatial 
SHS refers to correlations among local sub-regio ns within one 
footpo int si de of the magneti c arcad e, whereas lEmslie et all 
(2003.1) and IBattaglia & Benzl (l2006h compared spectra inte- 
grated over individual footpoint sources. We suspect that al- 
most all local regions in this event were subject to a common 
acceleration and transport process to share the same F-S re- 
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lation. For this reason we do not believe that our result is in 
conflict with the previous results. 

An entirely new result in this study is the spatial anticorrela- 
tion between HXR spectral index and electric field in the RCS. 
This property was found because we could measure both the 
spectral index and the electric field as functions of position 
within an extended HXR ribbon. As a comparison, we note 
that the numerical simulation for the direct electric field ac- 
celeration of electrons predicted electron energy distribution 
with power law index S = 4.0, 2.6, 1 .5 for el ectric field strength 
E = 0.1, 1.0, 10 V cm" 1 , respectively (IWood & Neukirchl 
l2005h . We are not sure how to convert this electron power 
law index S into the observed photon spectral index 7, be- 
cause further assumptions need to be made on the nature of 
the numerical solutions and the radiation. Within the scope 
of this paper, we note that there is a qualitative agreement 
betw een our observational r esult (Fig. [4]) and the model re- 
sult (IWood & Neukirchl l2005h . in that both indicate a hard- 
ening of electron energy distribution with increasing electric 
field strength. It appears that the present result suggests the 
dominance of direct electric field acceleration of flare elec- 
trons. It how ever does not exclude other possibilities. As 
discussed by Huds on & Farnfld (|2002|) . an explicit theoretical 
demonstration of the SHS behavior under the framework of 
the thick-target model can be ma de with the s tochastic accel- 
eration mechanism discussed by BenJ (|1977|) . Alternatively 
it is still possible that the efficient electron acceleration is 
confined in the regions of strong HXR flux and the weaker 
HXR flux regions nearby result from some propagation ef- 
fect. A more systematic survey of the relationship between 
these two physical quantities, both temporally and spatially, 
will be needed in order to ascertain such association under 
the context of electron acceleration mechanism. 
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